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Entangled light from driven dissipative microcavities 

D. PageQ and H. Feliske 

Institut fur Physik, Ernst-Moritz-Arndt-Universitdt Greifswald, 17487 Greifswald, Germany 

We study the generation of entangled light in planar semiconductor microcavities. The focus 
is on a particular pump configuration where the dissipative internal polariton dynamics leads to 
the emission of entangled light in a W-state. Our study is based on the nonlinear equations of 
motion for the polariton operators derived within the dynamics-controlled truncation formalism. 
They include the losses through the cavity mirrors, the interaction with lattice vibrations, and the 
external laser driving in a Langevin approach. We find that the generated entanglement is robust 
against decoherence under realistic experimental conditions. Our results show that pair correlations 
in solid-state devices can be used to stabilize the nonlocal properties of the emitted radiation. 

PACS numbers: 03.67.Bg, 42.50.Dv, 71.36.+C 


I. INTRODUCTION 

Quantum entanglement is known as the essential re¬ 
source for quantum information processing EH3I- It is 
defined as a nonlocal correlation that cannot be inter¬ 
preted in terms of classical joint probabilities gHH]. The 
most fundamental examples of nonequivalent forms are 
GHZ and W states jZj, [8]. The identification and quan¬ 
tification of entanglement is commonly based on entan¬ 
glement witnesses [shu]. The generation and control of 
entangled states still is a challenging task for quantum 
computation. In the optical domain, the implementa¬ 
tion of quantum algorithms relies on the availability of 
efficient sources for entangled photons. 

The generation of entangled photons is usually based 
on parametric down conversion in nonlinear crystals 
[H1IIS or biexciton decay in quantum dots Hang- Op¬ 
tically excited semiconductor microcavities [TT)H55] are 
alternative candidates for the efficient generation of en¬ 
tangled light on the /.xm scale [231428] . The exciting laser 
field with frequency near the fundamental band gap co¬ 
herently generates electron-hole pairs (excitons). The dy¬ 
namical evolution of excitons is governed by the Coulomb 
interaction, and the efficient coupling to the cavity pho¬ 
tons leads to mixed exciton-photon modes—so-called po- 
laritons PM5T] . The external laser can be tuned to stim¬ 
ulate parametric scattering processes between polaritons 
which may cause entanglement 1321 1331 . A moving polari¬ 
ton induces an electric polarization as a source of light 
that carries the initial (internal) polariton entanglement 
[22] . Then, depending on the explicit pump configu¬ 
ration, branch or frequency entanglement [231 l32| . po¬ 
larization entanglement [28] , multipartite entanglement 
[55] , and hyper-entangled photon pairs in multiple cou¬ 
pled microcavities [27] can be generated. 

In this work, we demonstrate that a semiconductor 
microcavity can be used to entangle light in a W state 
configuration. Beyond that, such a setup allows to an¬ 
alyze the internal polariton entanglement properties in 
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the presence of dissipation. Specifically, we consider a 
microcavity that is either continuously driven or excited 
by Gaussian pump pulses. In previous work [55], we in¬ 
troduced the specific pump arrangement for the creation 
of entangled light in a W state and used multipartite en¬ 
tanglement witnesses to verify the nonlocal correlations. 
Here, the inclusion of decoherence, induced by the losses 
through the cavity mirrors, and the coupling to lattice vi¬ 
brations, within the dynamics-controlled truncation for¬ 
malism, allows us to study the emitted light under realis¬ 
tic experimental conditions. Most notably, we show that 
the entanglement of the generated light is robust against 
dephasing. 

We proceed as follows. In Sec.[TT]we briefly recapitulate 
the equations of motion obtained within the dynamics- 
controlled truncation formalism and review the explicit 
pump configuration. The tomographic reconstruction of 
the state of the emitted radiation is performed in Sec. 
|III[ including the analytical solution in the limit of con¬ 
tinuous pumping in Sec. |III A| and the numerical solution 
for Gaussian pump pulses in Sec. mm Further details 
for the derivation of the analytical result can be found in 
App. 0 We finally conclude in Sec. [IV] 


II. THEORETICAL DESCRIPTION OF 
PARAMETRIC EMISSION 

The theoretical description of the dynamical processes 
in semiconductor microcavities is frequently based on 
an explicit bosonization of the whole system Hamilto¬ 
nian [22] ;24] [34]. Alternatively, one can derive equa¬ 
tions of motion for generalized Hubbard (transition) op¬ 
erators and truncate these equations at a certain order 
of the external field. This approach is called dynamics- 
controlled truncation scheme [351438] ; it can naturally be 
used to evaluate correlation functions that allow for a to¬ 
mographic reconstruction of the state of the emitted light 
modes [26] ;28j. It is thus well suited to study the gen¬ 
eration of entangled light in semiconductor microcavities 
under realistic experimental conditions. 
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A. Nonlinear system dynamics 


We begin with the equations of motion for the semi¬ 
conductor exciton and cavity photon operators that are 
derived within the dynamics-controlled truncation for¬ 
malism [371 EE], 

^a k = -iw k a k + if l R b k , (la) 

^6 k =-iw k 6 k + if2 fl a k -ii? k 1 ', (lb) 
where Rg L = R^ at + R:^, 

R i" = ^ E ‘L, +fa -A,»fe. (2»> 

ki,k 2 

K x = v xx J2 b L+k 2 -kK b k 2 ■ ( 2 b) 

ki,k 2 

In these equations, a k (6 k ) annihilates a cavity photon 
(semiconductor exciton) with in-plane wave vector k and 
energy (w k ). Qr is the dipole coupling strength 

between excitons and photons—the so-called Rabi fre¬ 
quency. In addition, n sa t denotes the exciton saturation 
density and V xx is the exciton-exciton coupling strength. 

The unitary Hopfield transformation to polaritons [29], 


PikA _ f 3fi k C lk \ /6 k \ 
V2k) \W 2k C‘2kJ \a k J 


(3) 


which is tuned to diagonalize the linear part of the equa¬ 
tions of motion ([l]), leads to the equations of motion in 
the polariton basis, 

j^Pik = -iwikPi k - i Rik , (4a) 

^;P2k = ^iw 2 kP2k - i^?2k L > (4b) 

with Rj\J‘ = XjkR^ L . Here, pj k annihilates a polariton 
with dispersion Wj k in the lower (j = 1) or upper (j = 2) 
branch. 


B. External driving and dissipation 

To include losses through the mirrors, the interac¬ 
tion with lattice vibrations and the external laser driv¬ 
ing, we couple the system dynamics to the environ¬ 
ment. As shown in Ref. m, combining the dynamics- 
controlled truncation scheme with the nonequilibrium 
quantum Langevin approach, the incoherent system dy¬ 
namics decouples from parametric scattering processes. 
In particular, we have to add the damping rates T^ 4 ’ and 
Langevin noise source operators T with proper statistics 
and moments to the equations of motion Q. This yields 

^Pik = -iwiKPifc - i Rik + Jp lk , (5a) 

T7P2k = -iw 2 kP2k - i^2k L + Xp 2k , (5b) 


where uij k = cc ? k — ir( ? 1 ( rf ' i /2. The operators T are charac¬ 
terized by vanishing expectation values, (Tp) = 0, where 
p = Pj k , and by the second order moments 

= 2(ZV(i)> 6{t - t') (6) 

with diffusion coefficients 

2 (D^(t)) = £)> - (AKM*) + vitMt)) ■ (7) 

In Eq. ([7]), the dot denotes the time-derivative following 
from Eqs. Q, i.e., without noise source operators T. 

Equations of motion for the expectation values 
—the so-called polariton photoluminescence— 
are given in Ref. [37] . They are derived in the framework 
of a second order Born-Markov approach. Important for 
us is the final result: Due to the decoupling of incoherent 
dynamics and parametric processes the diffusion coeffi¬ 
cients in Eq. © can be used as input when we calculate 
multi-time correlation functions of polariton operators. 
We stress that the damping rates follow from this 

treatment too. 


C. Explicit pump scenario 


Let us now consider the experimental setup intro¬ 
duced in Ref. [53], where four pump lasers drive the 
lower polariton branch at wave vectors k p i — ( kp,k p ), 
k p2 — ( A'rj 1 hp ), k p 3 — ( fcp, ^p), and kp 4 — (A: p , ^-p) 

(see Fig. [l|. The incident angles of all pumps are below 
the magic angle [hig&i such that single-pump scattering 
processes (signal at k and idler at 2k p „ — k) are negligible. 
The multi-pump parametric processes (signal at k and 
idler at k p „ + k pm — k with n ^ m) share a common idler 
mode at k^ = (0, 0). The four corresponding signal modes 
at k s i = (0,2fcp), k s2 = (-2fc p ,0), k s3 = (0,-2 k p ), and 
k s4 = (2fc p ,0) have been shown to be entangled | 33] . 

To obtain the equations of motion for the signal and 
idler modes, we introduce a simplified notation. In par¬ 
ticular, because all scattering processes are within the 
lower polariton branch, we omit the branch index. In 
addition, we introduce Y x = Y\ k;r for every quantity 
Y = P,u,T( tot \Z,R NL ,X,C and define 7x = T { x ot) / 2 
for x = i,s 1,..., s4,pl,... ,p4. Because of the particu¬ 
lar pump-signal-idler configuration, we have Lo sn = w s , 


7sn — 7 S 5 X sn — A s , Csn — ^pn 


and C pn = C p for n = 1, 


5 Xp n — Xp , 


,4. 


Assuming classical pump fields {pjk pn ) = Pn £ C, 
which imply a coherent driving, and identical pumps 
V n = V, we retain only terms containing the semiclas- 
sical pump amplitude V twice. Introducing the vectors 
P = (jPuph, ■ ■ ■ ,p\l) T and F = (F Pi ,F p ^,... ,F p ^) T , 
the equation of motion for the signal and idler modes 
takes the form: 


^P(t)=M(t)P(i)+F(t), 


( 8 ) 
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FIG. 1. Sketch of the considered pump-signal-idler config¬ 
uration in momentum space. The circles depict the possi¬ 
ble energy- and momentum-conserving scattering processes 
within the lower branch, where two pumped polaritons scat¬ 
ter into pairs of signal/idler polaritons. Specifically, mixed- 
pump scattering processes of oppositely arranged (neighbor¬ 
ing) pumps with |k p „ + k pm | = 0 (|k pn + k pm | = 2 k p ) con¬ 
tribute to the circle(s) with radius y/2k p (fc p ). 


where 


/ -i Wi 

-i 9sV 2 

-i g s V 2 

-i 9sV 2 

-i 9sV 2 \ 

i 9s(r*y 

• 5k 

1W S 

0 

0 

0 

i 9s(T*) 2 

0 

• ' * % 

1W S 

0 

0 

i 9s(V*) 2 

0 

0 

• 5k 

0 

Vi 9s{V*) 2 

0 

0 

0 

iw* / 


(9) 


and 

g s = 2X,X s X p (C p + V xx X p ) . (10) 

\ ^sat J 

We note that the matrix M depends on time solely, be¬ 
cause the pump amplitude V is time-dependent. 

D. State of emitted field 

Defining the matrix G of Green functions as the solu¬ 
tion of the homogeneous equation, 

^G(f,f') = M(f)G(M') (11) 

with initial condition G (t,t) = I (I is the 5x5 identity 
matrix), the solution of the inhomogeneous Eq. ([8]) is 

P(t) = G(t,0)P(0)+ f G(t, r)F(r) dr . (12) 

Jo 

It allows for the calculation of multi-time correlation 
functions. 

As a basis for the tomographic reconstruction of the 
measured signal/idler photon density matrix we choose 
the four states |lj,l sn ) (n = 1,...,4), where 11^) de¬ 
notes the state of a photon in channel x = i, si,..., s4. 
This choice can experimentally be realized by postselec¬ 
tion of events, where a click in the idler detector oc¬ 
curs, which takes out the vacuum component. Then the 
matrix elements of the measured photon density matrix 
Pi,sm;i,sn = (U, ^sm\p\U, Gn) are given by 


Pi, 


sm\i,sn 


h. 

h. 


’T d JT d 


’T d JT d 


(Pi(tl )p\ m (*2 )Psn (f 2 )Pi {t 1 ) ) dt i dt 2 

(pl(tl)Pi(.h))(plm(.h)P S n(t 2 )) + (pl(tl)pL(*2))(p S n(t2)Pi(il))l dt 2 , 


( 13 ) 


where the second line follows from a Wick factorization. A. Analytical modeling 

In this equation J\f is a normalization constant and Td is 
the detector window. 


III. RESULTS 

The equations from the last section allow us to study 
the tomographic reconstruction of the state of the emit¬ 
ted signal/idler fields in different situations. If the semi¬ 
conductor microcavity is continuously pumped, the equa¬ 
tions of motion can be solved analytically through trans¬ 
formation into the pump rotating frame. For Gaus¬ 
sian pump pulses—the usual experimental situation—the 
equations have to be solved numerically. 


To obtain analytical results for the stationary state 
in the long-time limit we assume a continuous pumping, 
i.e., V = 'Pe~ 1Upt with Pel. We define A = g s V 2 
for abbreviation and perform a transformation into the 
pump rotating frame: 


= 7}! e 

r sn rsn 1 


F t 

Pm. 


= F t 

Pin 


2i uj v t 


( 14 ) 
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Defining 


T(i) = 


1 

0 

0 

0 

Vo 


0 

a — 2iujpt 

0 

0 

0 


0 

0 

0 

r,—2i LOpt 

0 


0 
0 
0 
0 

-2i ujpt J 


(15) 


and 


P(f)=T(f)P(f), F(f) = T(f)F(f) (16) 
brings the equation of motion ® to the form 


A P (t) = MP(t) + F(t), 
with the time-independent matrix 


(17) 







/o 0 
0 1 

0 

0 

0 

0 

°\ 

0 

M = T _ 1 (t)M(i)T(f) - 

2i uj p 

0 0 

1 

0 

0 






0 0 

0 

1 

0 






0 

0 

0 

0 

1/ 


(~T 

-iA 

-iA 

-iA 

- iA v 





iA 

-Is 

0 

0 

0 




= 

iA 

0 

-Is 

0 

0 

— 

iw 

<1. 


iA 

0 

0 

-Is 

0 





\ iA 

0 

0 

0 

-Is) 





(18) 


According to G(t,t') = T 1 (t)G(t,t')T(t'), the Green 
functions become a matrix exponential 

G(f, t') = exp{M(f — t 1 )} = G (t — t'). (19) 

In order to calculate the populations and corre¬ 
lators in the long-time limit needed for the tomo¬ 
graphic state reconstruction, we assume the condition 
7 s 7 i > 4A 2 to be fulfilled. This guarantees that all 
eigenvalues of M have negative real parts, i.e., the 
corresponding Green functions converge in the long¬ 
time limit. In addition, we assume a uniform noise 
background N b , characterized by {X p {t)F p (t’)) = 

<■^4 (*)^t (*')> = °> ^ (t)-F„ B (0> = N b^ y S(t - t'), 
and ((f) Jpt (t')) = (N b + l)T x S XiV 5(t - 1') with x,y = 
i, si,..., s4. As shown in App. [Aj the tomographic re¬ 
construction is 


P = 



/I 0 0 0\ 
0 10 0 
0 0 10 
Vo o o i/ 


( 20 ) 


with X £ [0,1]. The state p in Eq. ( |2()[ ) is a mixture of 
a pure (fully entangled) W state and a (not entangled) 
identity state, where the parameter X is the weight of 
the fully entangled W state in the mixture. For X = 1, 
p is fully entangled. Contrariwise, p is fully separable 
for X = 0. Clearly the state p is entangled for any finite 
X > 0. In this sense, X can be taken as an entanglement 



measure, which quantifies the violation of a correspond¬ 
ing Bell inequality. 

Figure [2] shows X as a function of A and N b . We note 
that the parameter A is proportional to the pump inten¬ 
sity, which, however, is limited by the stationarity condi¬ 
tion 4A 2 < 7 i 7 s . Obviously, the fully entangled pure W 
state is obtained for vanishing noise background. This 
result is in accordance with the discussion in our previ¬ 
ous article [33], where losses through the cavity mirrors 
and the coupling to lattice vibrations are neglected. In¬ 
terestingly, even for a finite noise background N b > 0 the 
pure W state can be generated if the pump power is high 
enough. Lowering the pump power at fixed N b leads to 
a decrease of entanglement. 


B. Numerical solution 


Numerically, we can also study the case of Gaussian 
pump pulses. In practice, we solve Eq. (Ill for Gaus¬ 


sian pump pulses, having an intensity maximum at 4ps 
and a width of 1 ps, and calculate the populations and 
correlators to do the tomography. Thereby we choose a 
reasonable detection window of Td = 120 ps, allowing for 
technically feasible experiments with standard photode¬ 
tectors. Again, the tomographic reconstruction results in 
a state of the form (201, i.e., it is fully characterized by 


a single parameter X. The results for the amount of en¬ 
tanglement quantified by A as a function of the pump in¬ 
tensity for different environment temperatures are shown 
in Fig. [3] Compared to the analytical solution in the last 
section, we have set the uniform noise background N b to 
0. Noise enters the equations through the pump-induced 
photoluminescence [40] that depends on the temperature 
of the reservoirs. The choice N b = 0 is the reason why 
the numerical results tend to 1 for vanishing pump in¬ 
tensity. Contrary to the long-time behavior for contin¬ 
uous pumps, an increase in the pump intensity leads to 
a decrease of the entanglement. This behavior is even 
more pronounced for higher reservoir temperatures. The 
reason for this is a temperature-dependent background, 
created by the pump-induced photoluminescence, on top 
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FIG. 3. Entanglement of the numerically reconstructed signal 
density matrix as a function of the pump intensity for different 
environment temperatures. 

of which parametric scattering, i.e., entanglement gener¬ 
ation, takes place. Increasing the temperature at fixed 
pump intensity leads to a higher background at a fixed 
number of parametric scattering processes and hence to 
a lower degree of entanglement. Nevertheless, the gener¬ 
ated entanglement is surprisingly robust (see the range 
of X in Fig. [3]), even in the full simulation which includes 
the losses through the cavity mirrors and the coupling to 
lattice vibrations. 


IV. CONCLUSIONS 

We have studied the generation of multipartite en¬ 
tangled light in semiconductor microcavities within the 


dynamics-controlled truncation scheme. Including the 
losses through the cavity mirrors and the coupling to 
lattice vibrations, this formalism allows for a decoupling 
of the incoherent system dynamics (pump-induced 
photoluminescence) from the parametric scattering 
processes as the source of entanglement. Calculating 
particular multi-time correlations functions, the state 
of the emitted signal/idler fields is obtained through 
tomographic reconstruction. The resulting multipartite 
entanglement between the four signal channels for 
both continuous pumping and Gaussian pump pulses is 
robust against decoherence under realistic experimental 
conditions. This observation shows that the emitted 
photons carry the initial polariton entanglement. Since 
polaritons are quasiparticles composed of cavity photons 
and semiconductor excitons, they can sustain pair 
correlations over long times and distances inside such 
solid-state devices. In this sense, the emitted photons 
serve as a probe of the internal entanglement properties. 
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Appendix A: Explicit solution for continuous pumping 


Here, we evaluate the stationary populations and correlations in the long-time limit. We start with the diagonal- 
ization of the matrix M from Eq. (18). The eigenvalues of M are 


Ai — A 2 — A 3 — — 7 s — iw.j , 


A 4/5 = ('ft + % ± V (7i - Is) 2 + 16A 2 ) - i Vi. 


(Al) 


To simplify the notation, we introduce Q = \J (7 * — q s ) 2 + 16A 2 , A = — 7 S — i 07 , and A± = A — A 4/5 = ( 7 ,; — 7 S ± f2)/2 
with A + > 0 and A_ < 0. With these definitions, the matrix V of eigenvectors is 


1 

•\/2DA + |A_| 


/ 0 0 0 
-v/DA+|A_| -\/DA + |A_| -\/DA + |A_| 
0 0 y/DA+|A_| 

0 v/DA + |A_| 0 

V v/DA+|A_| 0 0 


- V / 2|MA + V^A+|A-I\ 
^2|A_|iA a/^A+IA 
v /2|A_|iA y^AY+iA 
v /2|A_|iA -\/2A+iA 

^2]X^iA i/AAjTiA ) 


V = 


(A2) 
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such that 


V -1 MV = AI- 


/O 0 0 0 0 \ 

'OOOO O' 
0 0 0 0 0 
0 0 0 A+ 0 

Vo 0 0 0 xj 


(A3) 


The matrix G of Green functions is given by 


G (i) = e At V 


with matrix elements 


n 

0 

0 

0 

0 

\ 



(Gi4t) 

G, 

s(t) 

Gi 

s(t) 

Gi 

s{t) 

Gi 

s (t)\ 

0 

1 

0 

0 

0 




Gl s (t) 

G s 

At) 

Gs 

s'(t) 

Gs 

At) 

G s 

s' (t) 

0 

0 

1 

0 

0 


V" 1 

= e At 

G*, s (t) 

G s 

s'{t) 

G s 

s(t) 

Gs 

At) 

G s 

s'{t) 

0 

0 

0 

3 —A+t 

0 




GlM 

G s 

s'(t) 

G s 

s'(t) 

Gs 

At) 

G s 

s'(t) 

Vo 

0 

0 

0 

e -A_ 



\G*iJt) 

G s 

s'{t) 

G s 

At) 

G s 

At) 

G s 

At) J 


(A4) 


= e = (7 < I*)* ^cosh - 7 ‘ ^ 7s sinh jt'j , 
Gi, s (t) = -2i^e - 5 ( 7 ,- 7 ,,:)t sinh^f, 


3 e ~i( 7 i~ 7 s)t 

G s , s (t) = - + --- 


/ , n 

I cosh —t + 


7i 


IT 1 si ” h 


(A5a) 

(A5b) 

(A5c) 

(A5d) 


G s , s '(t ) — G St s(t) — 1. 

Convergence of these functions requires 7 i +y s > fl, i.e., 7 i 7 s > 4A 2 . 

Introduction of the uniform noise background N b allows for the evaluation of the idler and signal populations in 
the long-time limit. This yields 


N°° = = J 2 e 27at {n^GI^t) + 4{N b + 1 ) 7 * |G»,. s(r)| 2 } dr 


Nb'Yi ls(li +7s) -4A 2 (Nb + 1)7 S 4A 2 


N Z = , lim (pl(t)Ps{t)) = 


7i +7s 

(■ N b + 1 H 


lil « - 4A 2 

A 2 


t—> OO 

and the correlators become 


7 i + 7s 7i7s - 4A 2 


= I™ (■ PHt)PAt )) = {Nb + 1)71 


A 2 


li + 7s lils - 4A 2 


3 N b 
4 

Nt 

4 


7* + 7s 7i7s - 4A 2 

Afc7s 7i (7i + 7s) - 4A 2 
4(7i + 7s) lils ~ 4A 2 

A&7s 7i(7i + 7s) - 4A 2 
4(7i + Is) lils - 4A 2 


7V-= lim (P}{i)P${t')) = 

t,t —> OO 

Finally, the tomographic reconstruction results in the state 


P = 


1 N%NZ' + | Nl 


4 N$N°° + | N% 



i(2 N b + l)7i7 s A 
( 7 ® + Is)(lils - 4A 2 ) 


N°°N™, + in; 
N°°N™ s + | Nl 



1 0 0 0\ 
0 10 0 
0 0 10 
0 0 0 1 ; 


such that the value of A' in Eq. (201 becomes 


X = 


ATOO AT OO 

1,1 s,s' 

N°°N~ ■ 


1^ 


00 12 
i,s I 


(A6) 

(A7) 

(AS) 

(A9) 

(A10) 

(All) 
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